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ABSTRACT 


An annular turbine cascade, designed for laser-Doppler 
velocimetry , was further modified to provide additional laser 
and pressure probe access. The purpose of the research was to 
devise laser anemometry measurement techniques in a confined 
annulus and improve the ability to compare numerical 
predictions with experimental results. Flowfield computations 
were completed using a viscous flow solver, with the numerical 
exit plane coincident with experimental measurement location. 
A data reduction program was written to transfer non- 
dimensional numerical output to the experimental coordinate 
system for comparison of total pressure ratio, Mach number 
and flow angle. Endwall two-dimensional laser measurements 
were obtained through a 1.0922 millimeter hole to a depth of 


approximately three millimeters. 
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I. INTRODUCTION 


The non-intrusive measurement of secondary flows in 
annular turbines can provide insight into the design and 
performance of turbomachinery. Turbine stage efficiency can 
be improved with more precise understanding of the flow 
mechanics and the ability to numerically predict the response 
to blade redesign without manufacturing and testing. 

Goldman and Seasholtz (Refs. 1-3} have conducted 
continuing research in the field of laser-Doppler velocimetry 
(LDV) measurement and numerical prediction in an annular 
turbine cascade. Their efforts have concentrated on obtaining 
comparisons ahead, in, and downstream of the blade passage. 
In order to obtain their results, an extensive optical system 
with correction for curved surfaces has been developed. To 
obtain two-dimensional velocity, they have measured one 
component ard then rotated the optics to resolve the second 
component. 

This report presents the results of laser-Doppler 
velocimetry techniques and three-hole pressure probe 
measurements in the wake of a high turning turbine nozzle. 
Only radial fiber-optic probe traverses were performed. 
Circumferential surveys were obtained by manually rotating the 
annular cascade blade row upstream of the measuring location. 


Pressure measurements were compared to numerical results 


obtained using a three-dimensional viscous computer program. 
Laser-velocimetry techniques using two-component fiber-optic 
equipment w2re developed for the small annulus’ cascade 
including two-dimensional measurements of the endwall flow 


through a small access hole in the outer casing. 


Tt. EXPERIMENTAL SETUP 


A. TEST FACILITY 

Airflow for the annular turbine cascade experiment was 
provided by a VA-312 Allis-Chalmers 12 stage axial-flow 
compressor housed in the Naval Postgraduate School 
Turbopropulsion Laboratory. The compressor was operated at 
12,000 rpm at various discharge pressures. This air supply 
was also used to generate up to Mach 0.795 through a 107.95 mm 
(4.25 in) diameter nozzle for calibration of a three-hole 
pressure probe. 

Upstream air pressure was metered through the use of two 
dump valves and one cutoff valve prior to a plenum chamber. 
Air from the plenum was directed to the test section through 
flow straigh‘tteners in a 254 millimeter (10 in) steel pipe as 
shown in Figure 1. Flow stagnation pressure was measured 
upstream at two locations. One combination probe was 
connected to a mercury manometer board located near the dump 
valve controls. This pressure-setting probe also provided a 
digital readout of flow stagnation temperature. The second 
probe, and the static ports in the air supply pipe, were 
connected to a Scanivalve for automated pressure acquisition. 

LDV seed material, atomized-glycerin particles, was 


introduced through an eight millimeter (5/16 in) diameter 


copper pipe as shown in Figure 1. Seed atomization was 
accomplished using laboratory compressed air with a commercial 
TSI, Inc., Six-Jet Atomizer. Laser traverse table buffet from 
the exiting flow was eliminated with a prefabricated aluminum 


cone (Figure 1). Flow characteristics were examined with the 


cone on and off, using pressure measurements and dye injection 


from the twc chord inner hub static ports. 





FLOW STRAIGHTENERS f 
\ 





> INNEN HUB STATIC PRESSURE TUBING GUIDE 


Figure 1. Supply Pipe Schematic 


B. ANNULAR TURBINE CASCADE 

1. GENERAL DESCRIPTION 
The Annular Turbine Cascade (ATC) was developed to 
provide a small-scale device for flow experiment and 
computation on a transonic turbine stator. Thomas [Ref. 4] 


developed the basic apparatus. Figure 2 shows a cutaway view 


of the construction of the ATC with a radial view of the blade 
passage just upstream of the test section. Figure 3 shows the 
blade geometry on a section perpendicular to the radial plane 
and Figure 4 shows an end view with the wake positioning 
capability. The inner hub radius was 98.81 mm (3.89 in) with 
constant blade geometry to the outer radius of 116.46 mm 
fios5 in). The stator row had 31 blades with a midspan 
Spacing of 21.77 mm (0.857 in). Reference 4 contains the 


original set of drawings used to manufacture the ATC. 


eT t 
le ’ y 


ih 








Figure 2. Annular Turbine Cascade Cutaway 


2. PRESSURE MEASUREMENT MODIFICATIONS 

The original ATC pressure measurement access was 
located two axial chords downstream of the blade trailing 
edges. Figure 5 details the old and new probe holes. The 
outer casing was modified to accommodate the half-chord access 
hole and is shown in conjunction with the adjustable holder. 
A new set of static pressure ports was added to the inner hub 
as shown. The ports were set 90 degrees apart and in axial 
alignment with the original holes. Tubing was connected to 
provide an average pressure tap for both hub axial locations. 
Appendix A contains a complete set of the drawings used for 


machining. 
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Figure 5. Pressure Measurement Changes 


3. LDV MODIFICATIONS 

LDV access was changed to minimize light reflection 
and improve off-axis options. The desire to ultimately obtain 
a third component necessitated an offset angle for a second 
fiber-optic probe. The design allowed for simple removal and 
replacement of the window during operation for cleaning, and 
radial access as shown in Figure 6. The position of the wake 
relative to the probe volume was controlled externally by 
rotating the blade row, inlet bellmouth and center body ahead 
of the measurement volume. An alignment blank was used to 
orient the een traverse axis to the ATC. The 1.0922 mm 
(0.043 in) diameter hole in the blank provided an endwall flow 


measurement capability as shown in Figure 7. 


2-D Laser Component 


1-D Laser Component 


' 

{ 
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| | 
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Figure 6. LDV Access Modifications 
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Figure 7. Endwall Measurement Technique 


C. PRESSURE PROBE 

Total pressure, Mach number and flow angle measurements 
were made with a three hole yaw probe [Ref 4]. All pressure 
taps were connected to a Scanivalve with a +344 kPa (+50 psi) 
differential transducer. Scanivalve control and data 
recording were accomplished using an HP-9000 computer system. 
Appendix B contains calibration curves anda listing of the 
Hewlett-Packard code used to acquire and output data. 
Pressure data from the ATC were acquired with the same 
hardware and software used in the probe calibration. 

Previously, yaw alignment was accomplished using a "U"'- 
tube water manometer. Difficulty encountered with hysteresis 
and the response of the water to yaw misalignment was 


corrected wixxh the use of a differential pressure transducer. 


D. LASER-DOPPLER VELOCIMETER 

The LDV system is shown schematically in Figure 8. The 
Lexel Model 95 four-Watt Argon-Ion LASER was connected to a 
TSI, Inc., Model 9201 ColorBurst Multicolor Beam Separator. 
The beam separator split the incoming light into two beams, 
with one beam passing through a Bragg cell for 40 Mhz 
frequency shifting. The two beams were further split into 
three polarized pairs: green (514.5 nm), blue (488 nm), and 
violet (476.5). The TSI, Inc., Model 9230 CoG mien 
Multicolor Receiver provided the frequency-shift signal to the 
ColorBurst for beam shifting. Each beam was directed through 
individual couplers for transmission to the laser probe by 
fiber-optic cable. The fiber-optic probes contained receiving 


optics to direct the return signals to the ColorLink. 





Figure 8. LDV Wiring Schematic 
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The ColorLink provided photomultiplier and frequency 
Batting functions. Control of shifting and photomultiplier 
voltage was accomplished remotely by computer. After 
conditioned by the cColorLink, the signals were sent to 
individual TSI, Inc., Model 1990C signal processors where 
valid Doppler signals were counted with reference to desired 
characteristics set by the user. An MI-990 interface 
monitored coincidence criteria (2-D measurements) and 
transferred valid counter information to computer memory. 

The fiber-optic probes were mounted to a Model 41583 
LINTECH Traverse Table. An Applied Motion Products System 
1618 Traverse Controller was used to control traverse movement 
manually and by computer. [Ref. 5] All remote control of LDV 
system equipment and data processing was completed through the 


menu driven environment of TSI’s FIND version 4.04. 


E. EXPERIMENTAL PROCEDURE 
1. PROBE CALIBRATION 

Probe calibration and investigation ue data 
acquisition options were completed together. The calibration 
was accomplished using a free-jet with a 107.9 mm (4.25 in) 
nozzle connected to the air supply (Figure 1) used for the 
annular cascade. Figure 9 shows the testing apparatus and 
pressures used to calibrate the cobra probe. 

The probe was aligned in yaw with a +68.9 kPa (+10 


psi) Stathan., SCO 277-10, differential pressure transducer. 


ig! 


Yaw-alignment response was an improvement over the previous 
manometer arrangement. The transducer bridge output was 
periodically zeroed by removing the probe from the flow. The 
yaw pressure measurement was tapped from one side of the 
transducer as shown in Figure 9. The Scanivalve system 
recorded the five pressures shown in Figure 9. Pre-operation 
tubing-leak checks were completed by trapping pressure to each 


Scanivalve port in turn. 








Yaw Differential 
Pressure Transducer 






Probe Tip 








Pa 
#2 P2/3 

#3 Free Jet Static 
#4 Free Jet Total 
Plenum Total 


Figu::e 9. Yaw Probe Calibration Equipment 


The probe was calibrated at Mach numbers from 0.20 to 
0.79 in steps of 0.1. Pitch sensitivity was investigated at 
twelve positions from minus ten degrees to plus eight degrees 
(Figure 9). Appendix B contains a description of the 


calibration variables and shows the sensitivity of the probe 
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to pitch angle. Sixth-order polynomials were fit to the 
calibration data. The zero-degree probe calibration was 
successfully validated (Appendix B) under different 
atmospheric conditions. 

2. PRESSURE PROBE MEASUREMENTS 

Wakes were surveyed by repositioning the blade passage 
upstream relative to the probe. The inner hub rotated with 
the blade ring as the wake was positioned for LDV or pressure 
measurement. The bolt shown in Figure 4 was used to position 
and hold the inner ring and stator. A total of 18 degrees 
circumferential variation, in steps of one degree, was 
achieved. 

Cobra-probe measurements were completed as described 
in the calibration section. Upstream static and stagnation, 
and inner-hub static pressures were used in conjunction with 
the center-port and yaw-port pressures of the probe. The 
probe provided total pressure and flow angle. Mach number was 
obtained from B (Appendix B). The desired hub-static to 
upstream-stagnation pressure ratio was set using’ the 
Scanivalve system and manual control of the air supply dump 
valves. The: probe was positioned radially and aligned with 
the flow before pressure data were acquired. Upstream 
stagnation temperature and wake flow angle were recorded 
manually for each measurement location. Flow angle was read 


from a vernier scale to an accuracy of plus or minus 0.2 


iS 


degrees at span locations greater than 25 percent. Span 
position was established with reference to 100 percent at the 
outer wall. 
3. LDV ALIGNMENT 

LDV alignment for endwall flow measurement was 
attempted as shown in Figure 10. The laser mounting bracket 
micrometer was set to ten millimeters for visual alignment 
around the access hole as shown. Obtaining a repeatable depth 
reference was not possible. This alignment technique did 
provide micrometer scale coverage for data recording and 
ensured the Doone passed cleanly through the hole. The length 
of the probe volume precluded the use of the outer surface of 


the plate as a reference point. 


1.0922 mm hole | 


Laser beams positianed at the edge of the access hole. 





Figure 10. Endwall LDV Alignment Technique 


14 


Alignment for automatic traverse of the wake passage 
through the window was accomplished using the alignment blank 
(Figure A8). The multi-component beams were positioned in 
each of the six alignment holes using manual control of the 
traverse table. The coordinates were read from the digital 
position indicator on the traverse controller for each hole. 
Axis conversion from the traverse to the ATC was completed 
using the pcsition data and the dimensions of the alignment 
blank (Appendix A). Traverse table movement could then be 
programmed using the automatic traverse control menu in TSI’s 
FIND software. 

4. ENDWALL MEASUREMENTS 

Endwall flow data were acquired without the use of 
automatic traverse control. The laser beams were aligned with 
the half-axial-chord access hole as described above, and then 
traversed manually using the probe-bracket micrometer. 
Adjustments to seeding position, counter gain and 
photomultiplier voltages were continually updated to maximize 
data rate and minimize the effect of noise. Seed-wand 
position was particularly important. Micrometer position 
information was married to the appropriate data files in post 
processing of FIND’s output ASCII .PRN files using Borland 


Corporation’s spreadsheet program Quattro Pro version 5.0. 
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IiIt. COMPUTATIONAL FLUID DYNAMICS 


A. GRID GENERATION 

Grid generation was completed using Turbomachinery C GRID, 
TCGRID, (Ref. 6). The code was developed to build 3-D "C" or 
"H' type grids for turbomachinery in the standard PLOT3D 
format (NASA Ames Research Center graphics’ software). 
Appendix C contains the namelist input file used to generate 
the stator g.zometry and a discussion of the output. The inlet 
bellmouth and center body as shown in Figure 4 were not 
modeled for this study. 

Careful manipulation of the namelist variables was 
required to »xbtain a smooth solution. The grid geometry was 
developed to improve resolution of the endwall boundary layer. 
The ability zo generate a fine grid was hampered severely by 
the length of wake desired in the model. Reduction of the 
wake-exit plane from one axial chord downstream of the blade 
trailing edges [Ref. 4] to one-half axial chord resulted ina 
dramatic increase in grid resolution. The final grid 
dimension was 150x31x65, up from 121x31x21. The blade surface 
received an 2dditional 29 points and the radial resolution was 
improved from 21 points to 65. The smaller computational wake 
length also reduced the interval between grid lines from the 


trailing edge of the blade. 
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B. COMPUTATIONAL SCHEME 

Flow solutions were completed using Rotor Viscous Code 3-D 
(eyC3D). RVC3D [Ref. 7] was designed to solve three- 
dimensional viscous flows in turbomachinery. The code solves 
thin-layer Navier-Stokes equations in Cartesian coordinates 
using an ex dlicit finite difference technique. Turbulence 
effects were modeled using the Cebeci-Smith option. The 
restart option of the code was utilized throughout the 
investigation to monitor convergence of the- solution. 
Appendix D contains the namelist input file used to obtain the 
final solution and a description of steps required to run the 
code on the Naval Postgraduate School (NPS) Cray computer. 
Output solution files (fort.3) were examined using the 
graphics software PLOT3D. Comparison with experiment required 
data reduction of the output file. Appendix E contains a 
synopsis of the output file structure with data reduction 
steps required for comparisons and a listing of the FORTRAN 
program (plajie.f) used to accomplish the manipulation. 

The fina’. solution selected for comparison was chosen for 
the minimum residual and the smoother transition of a solution 
across the wake cut. An area of unsteady flow became 
prominent near the hub and on further analysis revealed an 
area of predicted reverse flow. Figure 11 contains a 
comparison between total pressure ratio results at the 50 


percent span location for 1200 and 1490 iteration solutions. 
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One blade wake is shown with the wake cut in the grid clearly 
evident in the plot. The closely-spaced points in the wake 
cut of the TCGRID geometry has introduced a mathematical 
boundary. The 1200 iteration solution exhibited the beginning 
of the reverse flow area near the hub, but contained serious 
discontinuities because of the solution scheme. The 1490 
iteration total-pressure-ratio plot demonstrates the smoothing 


that took place as the solution progressed. 


CFD Solution Comparison 
Total Pressure Ratio (50 % Span) 


ae 34 -32 -30 
CFD Wake Position - theta (degrees) 


Iterations —— 14990 «© 100 | 


Figure 11. CFD Solution Comparison, 50% Span 
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IV. RESULTS AND DISCUSSION 


A. PRESSURE PROBE 

Annular Turbine Cascade pressure measurements were 
obtained at an average Reynolds number of 650000 based on 
axial chord. The pressure ratio, inlet total to exit static 
at the hub, was set at 0.68 to enable a proper comparison to 
be made with the CFD solution. Pressure and flow angle 
measurements were repeatable, with Scanivalve pressure- 
acquisition standard deviations of less than 0.20 percent. 
However, data obtained inside 25 percent span were at low 
Gynamic pressure and too difficult to resolve with the 
pressure probe, and are therefore not reported. Flow 
unsteadiness increased as the probe was immersed further in 
from the encwall. The direction of approach to zero yaw in 
the less-than-25-percent span region affected the flow angle 
measurement by as much as 0.5 degrees. 

Pressure and flow-angle data were taken at 18 wake 
positions for 25, 50, 75 and 90 percent span locations on the 
half-axial-chord wake measurement plane. Pressure data 
repeatability was demonstrated along the 50 percent span 
station for all wake positions, and radially at the two-chord 
zero-degree wake position. Figure 12 shows the total- 


pressure-ratio data obtained for the 50 percent span location. 
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The data show that there was a negligible influence of the 
conical diffuser. Figure 13 shows the repeatability obtained 
in radial surveys for a constant wake position at the two- 
chord measurement location. The repeatability shown in 
Figures 12 and 13 clearly validated the ability to reset the 
desired operating pressure ratio for the experiment. Flow 
angle meaSurementS were repeated at the 90 percent span 
location, and radially at the two-chord wake measurement 
location as shown in Figures 14 and 15. In all cases, the 
repeatability of the probe pressure meaSurements was 
excellent, with one largest discrepancy less than 0.5 percent 


in total pressure. 
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Figure 12. Pt/PO Repeatability, 50% Span 
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Figure 13. Pt/PO Repeatability, Radial 
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Figure 14. Angle Repeatability, 90% Span 
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Figure 15. Angle Repeatability, Radial 
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B. LASER MEASUREMENTS 

Two-dimensional LDV measurements were obtained through a 
1.0922 mm (6.043 in) hole to a depth of approximately three 
millimeters. Initial measurements were taken with one 
component (green) of the two-component probe. Figures 16 and 
17 show the velocity profile and turbulence intensity results 
from the first attempt at endwall measurements. The ‘depth’ 
shown in the graphs is an estimated distance from the endwall 
based on the rapid loss of data rate observed at the depth 
shown as 0.1 mm from the endwall. The data map a boundary 
layer profile to approximately 0.6 millimeters. The high 
pressure ratio (0.98) provided near perfect seed conditions, 
and the casing vibrations were low. The data shown for V-z do 
not necessarily coincide with the exact locations of the data 
for V-0 because the fiber-optic probe was rotated 90 degrees 
and alignment could not be maintained. 

The turbulence intensity results are with respect to the 
maximum velocity at the furthest measurement location into the 
annulus. The significantly higher turbulence intensity in the 
V-0 direction is probably due to the steep gradient of the v-6 


mean-flow profile in the radial direction. 
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Figure 16. 
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Figure 17. Single Component LDV Turbulence 
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Two conponent LDV measurements were obtained at 
approximately 15 percent higher flow velocity. Figure 18 
depicts a test section ’sweep’ of near-endwall velocity 
measurements at a pressure ratio of 0.974. Depth reference 
(in mms with respect to an unknown datum on a vernier scale) 
is shown. The exact distance from the wall was not 
obtainable. Flow periodicity is suggested and decreasing 
velocities toward the outer casing show the endwall boundary 
layer. The }moundary layer is distorted due to corner vortices 
and secondary flows that originate in the blade row. Flow 
angle and turbulence intensity plots are included in Figures 
19 through 21. Once again the turbulence intensity in the V-6 
direction is higher than that in the V-z direction. Figures 
18 through 21] represent automatic editing of the histograms by 
the FIND software (user selection). Data outside two standard 
deviations were discarded. 

Increased flow velocity introduced higher vibration levels 
and temperature increases as the test progressed. Hole 
alignment was periodically verified and, occasionally, small 


adjustments (less than 0.01 mm) were required. 
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Figure 18. Two Component LDV Endwall Velocity fae 
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Figure 19. Component Flow Angle 
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Two Component LDV Turbulence, 


Figure 20. 
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Figure 21. Two Component LDV Turbulence, V-z 
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C. NUMERICAL RESULTS 

The CFD solution was run with the same pressure ratio 
(0.68) used in the experiment. Graphical output of the 
numerical results was useful in obtaining a qualitative 
understanding of the flow characteristics. Figure 22 depicts 
the predicted pressure ratio for the blade surface at the hub, 
mid-span and tip. The suction surface curves suggest that the 
position of the blade passage throat was at 0.8 of axial 
choma. Figure 23 shows, with colored contours, pressure 
coefficient information for the blade surface pressure 
distribution. The Mach number contours shown in Figure 24 
provided an initial indication of where flow characteristics 
were well established. Boundary layers appear as expected, 
thin near the endwall and thick near the hub. Wake curvature 
in the spanwise direction is also evident in the Mach contour 
plot. Difficulty in obtaining repeatable flow angles in the 
experiment below 25 percent span may have been due to the 
unsteady fi.ow separation which 1S suggested by the 
computationa.. results for the region near the hub wall. 
However, the reverse flow predicted by the code in this region 


could not be verified in the experiment. 
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Figure 23. PLOT3D Blade Pressure Coefficient Plot 
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Figure 24. PLOT3D Exit Plane Mach Contour Plot 
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D. COMPARISON 

Three comparisons were made between numerical and 
experimental results. To show the flow periodicity , the CFD 
solutions were repeated over two blade passages as shown by 
the solid line in the following Figures 25 through 36. 
Comparisons between experiment and numerical (Appendix E) 
results are presented for total pressure ratio, Mach number 
and flow angie. Figures 25 through 28 show the total pressure 
ratio comparison. It is noted that uncertainty in the probe 
measurement of total pressure is small. The probe calibration 
uncertainty enters only when comparing Mach number. Wake- 
profile phase coincidence was established using the 50 percent 
Span Mach number data. The phase was kept constant for the 
circumferential comparisons at all other spanwise locations. 
Mach number (Appendix B) comparisons are contained in Figures 
29 through ~?2, and Figures 33 through 36 contain the flow 
angle data. 

CFD total pressure ratio predictions were within 10 
percent of the experimental values at 25 percent span. The 
difference decreased to below one percent at 90 percent span. 
Mach number comparisons with measurements were within 15 
percent at 25 percent span and decreased to three percent at 
90 percent span. Flow angle comparisons were within five 
degrees at 90 percent span with improvement to within two 


degrees at 25 percent span. The reason for closer agreement 


SZ 


in total pressure ratio and Mach number toward the tip and 
improving agreement in flow angle toward the hub is 


unexplained. 
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Figure 25. Total Pressure Ratio Comparison, 253 | 
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Figure 26. Total Pressure Ratio Comparison, 50% 
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Figure 27. 


Total Pressure Ratio Comparison, 75% 
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Figure 28. Total Pressure Ratio Comparison, 90% 





35 


Mach Number 
Comparison (25 % Span) 


Bo 
: 
a 
¥o 


1 
Normalized Wake Position 





Mach Number 
Comparison (50 % Span) 





Bo. 
3 

= 

Z 

3 

C 0.5 1 15 
Normalized Wake Position 
Figure 30. Mach Number Comparison, 50% a 


36 





Mach Number 
Comparison (75 % Span) 


weer er eet e® we eeeene2nee eee ee ew ew wee wee ewe ee ewe we ww eww wow we we ew ew we ewww Oe Be we wm wwe wm wn wwe ewww eee ee wee wee 
weerreer ere eeerenzeeeeeeeeee eee eee Fe me ewe wee ew ee ew ewww ew wet we ew ew ewww wee wnwe2eee2e eee e2ee2eeeeceeeeweecne 
werereeeeeeeee ene ew eeeer ee eee teen er ewe emer wee een eee eee ee eee eeewneneneneeece2nenn eee eee eceen eee 


eeeeeenere2eee2ceweeeeceereee2e ec ee eee ee ee we wee www www wwe wwe ew we eww ew eww e es wee ee ee ewe eee wow eee wees 


Normalized Wake Position 





Figure 31. Mach Number Comparison, 75% 
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Figure 32. Mach Number Comparison, 90% 
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Wake Flow Angle 
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Figure 33. “Flow Angle Comparison, 25% 
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Figure 34. Flow Angle Comparison, 50% 
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Wake Flow Angle 
Comparison (75 % Span) 
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Figure 35. Flow Angle Comparison, 75% 
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V. CONCLUSIONS AND RECOMMENDATIONS 


Experimental measurements and computational output 
locations were made to coincide. Repositioning the 
measurement volume to one-half axial chord downstream from the 
blade trailing edges improved the computational grid 
resolution. Flow solutions with the revised grid more 
realistically resolved endwall boundary layer characteristics 
so that future endwall LDV measurements will have a 
corresponding computed solution with equivalent resolution. 

Pressure probe surveys will continue to be important in 
comparing experiment with CFD prediction. The cobra probe 
calibration obtained was reliable and repeatable. An 
extensive set of baseline data has been accumulated for future 
use. The Scanivalve data acquisition system demonstrated the 
ability to obtain pressure readings with a standard deviation 
of less than 0.2 percent. The HP-9000 output was manually 
entered into a spreadsheet program for post processing. The 
time required to manually enter the output was overshadowed by 
the spreadsheet data manipulation flexibility. 

The most definitive comparison between CFD and probe 
measurements was in the total pressure ratio. Probe 
calibration error was not present in this comparison. The 
total pressire ratio disagreed by ten percent at the 25 


percent spar location and less than one percent at the 90 
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percent location. In all cases, the computer solution 
predicted a lower ratio. The difference may be attributed to 
computational modeling parameters or bias introduced in the 
experimental pressure ratio setting. Upstream total pressure 
was measured at only one station and this value may not have 
been a proper reflection of the average total pressure. Mach 
number compzrison characteristics followed the trend of the 
total pressu,sr‘e ratio as expected. Flow angle comparison with 
CFD was best at 25 percent span with a maximum difference of 
5° at the 9C percent span location. Of all the measurement 
locations, the 90 percent span location was the best defined 
and flow angle measurements made there were repeatable within 
0.2 degrees. 

Endwall LDV through a small (1.0922 mm) caSing access hole 
has been demonstrated. The ability to resolve endwall 


boundary layers and eventually tip flow losses appears to be 


promising. Precise LDV probe-volume positioning techniques 
were not obtained. The maximum pressure ratio attained for 
endwall LDV was 0.974. CFD and pressure probe data were 


accumulated at an operating pressure ratio of 0.68. Low data 
rate, vibration of the endwall hole, and thermal expansion 
problems will have to be dealt with as the operating pressure 
ratio is decreased to 0.68. 

Future efforts with the computational fluid dynamics 
should include modeling the inlet bellmouth and center body. 


The parameters selected for use in RVC3D should also be 
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explored to obtain a qualitative feel for cause and effect on 
the solution with particular attention to the unsteady region 
of the exit flow. 

Pressure data acquisition can be continued to obtain a 
more complete map of the wake~- region. Pressure 
instrumentation of a blade and passage would provide an 
additional source of CFD code validation data. Data reduction 
and presentation would be improved with a change in operating 
system. The current Scanivalve system with the HP-9000 does 
not allow file transfer to other computers for post 
processing. HP-9000 limitations restrict output flexibility. 

The annular turbine cascade now has off-angle LDV options 
that have not been explored. Three-dimensional radial LDV 
surveys of the wake passage should be attempted for comparison 
with pressure measurements already obtained. Initial 
indications show that the LDV probe-volume energy may not be 
able to provide return signals of sufficient intensity as the 
flow velocity increases. The addition of corrective optics to 
the fiber-optic probe may be required. 

The endwall technique developed in this study did not 
provide pressure equalization across the access hole. The 
pressure across the access hole should eventually be equalized 
to eliminate flow through the hole and therefore any influence 
the opening might have on the flow structure [Ref. 9]. The 
hole acted as a natural optical noise filter and further 


investigation of hole diameter influence on LDV performance is 
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warranted. In order to truly understand the hole influence 
and measurement reliability, the glycerin seed particle 
distribution and a study of their momentum characteristics 
should be undertaken. With these refinements, access to the 
flow througrk a hole in the casing is feasible, without 


changing the endwall flow characteristics. 
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Adjustable Pressure Probe Holder 


Figure Al. 
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Figure A2. 
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Outer Ring Modifications 


Figure A3. 
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Outer Ring Modification, Section A 


Figure A4. 
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Outer Ring Modification, Section B 


Figure A5. 
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Outer Ring Modification, Section C 


Figure A6. 
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Figure A7. 
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APPENDIX B. COBRA PROBE CALIBRATION 


The following synopsis was extracted from Appendix A of 
Reference 8. Mach number or dimensionless velocity can be 


related to a non-dimensional pressure coefficient, £6, as 


follows: 


Total velocity: Vicon PAGER 
Dimensionless velocity: A= —— 


Dimensionless velocity and Mach number are related as follows: 


CA Cole 
oe 
Jf =1-x? 
L. 
ya WYRT 
pit 
so that 
Y-iy 
x?=__4 
te ye 
2 
or 
fee 2X? 
yas (iL =X") 


BZ 


A pressure coefficient for the probe aligned in yaw can be 
written as 


Sof PSPo Ze) 2 Fa 


where p, and bp, are probe sensor measurements. Using p, and p,, 





respectively, it can be shown that 


g = 22 P23 _ (mM, Pitch) 
P, 


in the range of velocities which are of interest. 

Figure Cl contains the calibration data without curve 
fitting and Figure C2 shows the validation data plotted 
against the sixth order and linear curve fits at zero pitch 


angle. The sixth order curve fit equation is 


X=(a0+alB+a2B?+a3B?+adB4+a5Bh°+a6B*°)x100000 


with the coefficients shown in Figure Bl. 
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Cobra Probe Calibration Coeff 


Figure Bl. 
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Cobra Probe Calibration Data Plot 


Figure B2. 
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Probe Validation 


Figure B3.) 
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Data acquisition was accomplished using an HP-9000 


computer system with the following program: 


10 !FILE NAME : COBRA ATC 


20 !DISK LABEL "/AUSTIN" 

SOs LY JEFF AUSTIN 

40 !LAST MODIFIED 02/03/94 BY LCDR JOE SPITZ 

50m 

60 ! THIS PROGRAM RECORDS AND REDUCES MEASURED PRESSURES 
70% FROM A Scanivalve CONNECTED TO A THREE HOLE PRESSURE 
80 ! PROBE. Scanivalve PORT CONNECTIONS ARE AS FOLLOWS: 
90 ! 

100 ! PROGRAM VARIABLES AND PORT CONNECTIONS... (#3 not used) 
110 ! V = DESIRED S/V 

120°! A = PRESENT S/V PORT 

140 ! C = CHANNEL 

150 ! VO = DVM READING 

160 ! S$ = SCANNER LISTEN CODE 

170 ! S = SCANNER # 

180 ! Scanivalve #5 SET TO ADVANCE ONLY 

190 ! Pamb = Ambient Pressure (psi) 

200 ! Patm = Atmospheric pressure...Port # 1 

210 ! Peal = Calibration pressure...Port # 2 

220 ! Ptot = Reference total pressure...Port # 4 
230 ! Phub = Inner hub static pressure...Port # 5 
240 ! Pstat = Upstream static pressure...Port #6 
250 ! Pl = Center hole on probe...Port # 7 

260 ! P23 = Yaw port pressure...Port # 8 

270 ! beta = (P1-P23)/P1 

271! prat = Phub/Ptot (pressure ratio) 

272 ! mach = mach(P1, Phub, Gamma) 

280 ! 


290 OPTION BASE 1 

300 Gamma=1.4 

310 CLEAR SCREEN 

311 PRINTER IS CRT 

320 DISP "PLEASE WAIT WHILE RESETTING Scanivalves" 

330 PRINT 

340 V=5 

350 A=1 

360 GOSUB Read 

370 VK KKKKKEKREKKKEKKEKKKEKKKEKKKEKEKRKEKEKREKEKKREEKREEKEREEKKREKKKKKKKKKKK 
380 PRINT "DATA ACQUISITION FOR THREE HOLE YAW PROBE" 

390 DIM Y(20) 

400 DIM Patn(10,10), Pcal(10,10) , Ptot(10,10) , Phub(10,10) , Pstat (10,10) 
410 DIM P1(i0,10), P23(10,10) 

420 ! 

430 INPUT "ENTER MONTH, DAY, YEAR (IE. 02,05,94)" ,Y(3), ¥(5), Y¥(7) 


Di/ 


431 
432 
440 
441 
450 
451 
460 
550 
560 
570 
571 
572 
580 
581 
582 
600 
610 
630 
640 
641 
650 
660 
661 
662 
670 
671 
680 
700 
710 
720 
730 
740 
750 
840 
850 


860 : 


870 
880 
890 
900 
910 
ot2 
913 
914 
915 
916 


919 
220 
921 
ee 


INPUT “ENTER RUN #: ", Y(9) 

PRINTER IS 711 

PRINT USING "K,DD, rE 7 DD eset Dp UDATE OF RUN:",Y(3) ,¥(5) ,¥(7) 
PRINT "DATA RUN ",Y(9) 

PRINT 

PRINTER IS CRI 
PHAKKKKKKKKKEKKEKKKEKEKKKKKEKEKKEKKEREKKKKKKEKKKKKKEKKKKKEKKKKKEKKKKEKKEKKKEKK 
PRINT 

PRINT "CONSTANT FLOW. WAKE: +9 TO -9 DEGREES. SPAN: 5 LOCATIONS 
PRINT 


PRINT " Zero DVM on Scanivalve #5, port #1." 
PRINT 

PRINT " Set 20 inHg for Scanivalve calibration." 

PRINT 

PRINT" Adjust span for DVM = 0.009823 volts on port #2." 
PRINT 


PRINT" ENSURE DVM ADJUSTED AND DESIRED PRESSURE RATIO SET..." 
FOR WAKE = 1 TO 10 

INPUT “ENTER AMBIENT PRESSURE IN inHg", Pinhg 

Pamb=Pinhg* (.4911541) 

PRINT 


PRINT "CHECK DESIRED SPAN AND FLOW ALIGNMENT SET..." 
PRINT 
PRINT ''****kPRESS CONTINUE WHEN READY TO TAKE DATA****#" 
PAUSE 
CLEAR SCREEN 
PRINT 
FOR SPAN = 1 TO 10 
V=5 
A=1 
GOSUB Read 
PHKKKKKKKKKKKKKKKEKKEKKKKEKKKEKEKERKEKKKEEKKEKKKKKKKEKKRKEKKKEKKKKKKKKK 
/ 
! KKEKKEKKEKEKKKKKKKKEKKEKKKKKEKKEKKEKEKKEKKKEKKEKEKKKKKKKEKEKEKKRKKKKKKRKKS 
!RECORD THE PRESSURE DATA FOR THIS WAKE AND SPAN 
| 


OUTPUT 722;"F1R1T1M0Z1 " ! SETS UP HP3456A DVM 
V=5 ! Scanivalve #5 DESIGNATED 
S=1 ! HP3495A SCANNER #1 DESIGNATED 


PRINVTER IS 711 

PRINT 

PRINT "Wake Data Group ", WAKE 

PRINT "Span Location Counter ", SPAN 

PRINT USING "/,5X,K,DDD,/,2X,K,15X,K,17X,K";" 
Scanivalve #",V,"Port", "Volts", "Psi" 

PRINT 

PRINTER IS CRT 


! REZD Scanivalve PORTS 
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923m) 
930 
940 
950 
960 
970 


980 

990 

1000 
1002 
1010 
1020 
1021 
1022 


1023 
1024 


1025 
1026 
1027 
1030 
1031 
1032 


1033 
1034 


1035 
1036 
1037 
1040 
1041 
1042 


1043 
1044 


1045 
1046 
1047 
1050 
1051 
1052 


1053 
1054 


1055 
1056 


FOR A=1 TO 8 
GOSUB Read 
WAIT 2.0 


OUTPUT 701;"C" 
OUTPUT 701 USING "DDD";V+9 


! CLEARS SCANNER #1 


PORT 


TRIGGER 722 
ENTER 722;V0O 


| dere He de He He He He He HK KH HK KK HK KH HK HK HK KK KKH HK HK KI IKI KKK KK KEKE KEE KEKE 


PRINTER IS 711 


Pressure in psi 


IF A=1 THEN 

Patm (WAKE, SPAN) =V0*1000 

PRINT USING "2X,DDD,7X,7DD.DDDDDDD, 8X, 7DD.DDDD"; 
P, VO, Patm (WAKE, SPAN) 

PRINTER IS CRT 

PRINT USING "2X,DDD,7X,7DD.DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Patm(WAKE, SPAN) 

PRINTER IS 711 

END IF 


IF A=2 THEN 

Pcal (WAKE, SPAN) =VO*1000 

PRINT USING "2X, DDD, 7X, 7DD.DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Pcal (WAKE, SPAN) 

PRINTER IS CRT 

PRINT USING "2X, DDD, 7X, 7DD.DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Pcal (WAKE, SPAN) 

PRINTER IS 711 

END IF 


IF A=4 THEN 
Ptot (WAKE, SPAN) =(VO*1000) —Patm(WAKE, SPAN) +Pambb 
PRINT USING "2X,DDD, 7X, 7DD. DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Ptot (WAKE, SPAN) 

PRINTER IS CRT 

PRINT USING "2X, DDD,7X, 7DD. DDDDDDD, 8X, 7DD. DDDD"'; 
P, VO, Ptot (WAKE, SPAN) 

PRINTER IS 711 

END IF 


IF A=5 THEN 

Phub (WAKE, SPAN) =(VO*1000) —Patm (WAKE, SPAN) +Pamb 

PRINT USING "2X, DDD, 7X, 7DD.DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Phub (WAKE, SPAN) 

PRINTER IS CRT 

PRINT USING "2X, DDD, 7X,7DD. DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Phub (WAKE, SPAN) 

PRINTER IS 711 

END IF 
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! SCANNER SET TO READ Scanivalve 


IF A=6 THEN 
Pstat (WAKE, SPAN) =(VO*1000) ~Patm (WAKE, SPAN) +Pamb 
PRINT USING "2X,DDD,7X, 7DD.DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Pstat (WAKE, SPAN) 

PRINTER IS CRT 

PRINT USING "2X, DDD,7X,7DD.DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, Pstat (WAKE, SPAN) 

PRINTER IS 711 

END IF 


IF A=7 THEN 

P1 (WAKE, SPAN) =(VO*1000) —Patm (WAKE, SPAN) +Pambb 

PRINT USING "2X,DDD,7X, 7DD. DDDDDDD, 8X, 7DD. DDDD" ; 
P, VO, P1(WAKE, SPAN) 

PRINTER IS CRT 

PRINT USING "2X,DDD,7X,7DD.DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, P1 (WAKE, SPAN) 

PRINTER IS 711 

END IF 


IF A=8 THEN 

P23 (WAKE, SPAN) =(VO*1000) -Patm (WAKE, SPAN) +Pamb 

PRINT USING "2X, DDD, 7X, 7DD. DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, P23 (WAKE, SPAN) 

PRINTER IS CRT 

PRINT USING "2X,DDD,7X,7DD. DDDDDDD, 8X, 7DD. DDDD"; 
P, VO, P23 (WAKE, SPAN) 

PRINTER IS 711 


prat=Phub (WAKE, SPAN) /Ptot (WAKE, SPAN) 
beta=(P1 (WAKE, SPAN) —P23 (WAKE, SPAN) ) /P1 (WAKE, SPAN) 
mach=SOQRT (5* (P1 (WAKE, SPAN) /Phub (WAKE, SPAN) ) * (1/3.5)-5) 


END IF 

NEXT A 
PRINT "prat = ", prat 
EPRiieebeta = ", beta 


PRINT "mach # =", mach 
PRINTER IS CRT 
PRINT 
PRINT "Set Span and align probe. ..CONTINUE WHEN READY..." 
PRINT 
PAUSE 
CLEAR SCREEN 
NEXT SPAN 
PRINT "Verify Wake & update Pamb...CONTINUE WHEN READY..." 
PRINT 
PAUSE 
CLEAR SCREEN 
NEXT WAKE 
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4 

PKK KKKKKKKKKKEEKEKKEND OF MAIN PROGRAMK%& RR KK kk tok tok 
] 

! SUBROUTINE TO POSITION AND READ S/V PORTS 

! Scanivalve #1 HAS INOP HOME FEATURE 
Read: OUTPUT 707 USING ’#,K";V 

PO=SPOLL (707) 

L=BINAND(PO,15) 

T1=ROTATE (PO, 4) 

M1=BINAND(T1,7) 

P=10*M1+L 

CLEAR 707 

IF P=A THEN Finish 

OUTPUT 701; "Cc" 

OUTPUT 701 USING "DDD";V-1 

OUTPUT 701; "Cc"! 

WAIT .1 

GOTO Read 

Finish: RETURN 

END 
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APPENDIX C. GRID GENERATION 


Grid generation was completed on IRIS workstations. 
TCGRID, a FORTRAN program, required the adjustment of 
parameter statements to dimension arrays for the model build. 
Parameter aditustment required patience. The author found the 
manual helpful in most cases, however, the most useful tool 
turned out to be a running summary of the namelist input for 
troubleshooting. The most difficulty was observed in 
obtaining a smooth solution on the suction side of the blade 
geometry. The following TCGRID namelist file was utilized to 
obtain the grid shown in Figures Cl and C2. Reference 4 
contains the original list of blade coordinates and is not 
included here. 

&naml im=150 jm=31 km=65 itl=20 icap=18 k2d=3 merid=0 &end 

&nam2 nle=16 nte=10 dsle=.018 dste=.003 dshub—.0004 dstip=.0004 

dswte=.001 dswex=.060 dsthr=1. dsmin=.0004 dsmax=.025 
dsra=.45 rcorm.098 &end 

&nam3 iterm=100 idbg=0 00000 0 O aabb=.5 ccdd=.45 &end 


&nam4 fee 0 1-5 Oe 5a ee as. 
moee— 93.095 3.895 3.895 4.585 4.585 4.585 &end 
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Figure Cl. Grid Geometry On Inner Hub Surface 
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Figure C2. 
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APPENDIX D. RVC3D CODE INPUT 


RVC3D runs were initially completed on IRIS workstations 
at NPS. Solutions were limited to less than 500 iterations to 
ensure residuals were decreasing. File transfer of the fort.1 
(grid) and the rvc3d.in (namelist) files to the Cray were 
completed only after initial testing was completed. The 
following namelist contains the commands used for the final 
set of iterations to the minimum residual attainable. It is 
a restart (iresti=1) from 1200 iterations with termination at 


1490 iterations (imax=290). 


‘TRANSONIC TURBINE Annular Cascade - (Restart 12-1490) ’ 

&n11 im=150 jm=31 km-65 itl=20 i11=67 &end 

&nl2 cl=5.5 aviscl=0.0 avisc2=0.0 avisc4=0.50 ivdt=1 
nstg=4 itmax=290 irs=1 epi=0.50 epj=0.60 epk=0.60 &end 

&nl3 ibcin=3 ibcex=3 isymt=0 ires=10 icrnt=50 
ilresti=1 iresto=1 ibcpw=0 iginm0O &ernd 

&nl4 emxx=0.13 emty=0.0 emrz=0.0 expt=0.0 prat=0.6800 
ga=1.4 om=0.000000 igeom=1 alex=-67.0 &end 

&nl5 ilt=3 tw=1.00 renr=6.651e6 prnr=.7 prtr=.9 
Vispwr=.666666 srtip=0.0 cmutme14. jedge=31 kedge=20 
lltin=2 dblh=0.0024 dbl1t=0.0089 &end 

&n16 101=1 102=165 oar=0. ixjb—0 njo=l1 nko=3 jo=1 
ko=5 11 16 &end 


Runs on the NPS cray were accomplished with the command 
>qsub myjob.ngs 
where the file myjob.ngs was as follows: 
myjob.ngs 
#QSUB -q reg -1T 43200 -1M 9Mw 


cd cfd 
rve3d <rve3d.ine Lve3a- our 
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The RVC3D executable, grid (fort.1), restart (fort.2), 
and the input mnamelist (rvc3d.in) were held in_ the 
subdirectory cfd. The final solution required 13.2 hours of 
Cray computer time. A plot of the residuals is shown in 


Figure Dl. 


RVC38D Residuals 


Normalized to Maximum 


RMS(1)/1.59E-05 


GOO 800 1000 
Iteration 





Figure Dl. RVC3D Residuals 
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APPENDIX E. RVC3D DATA REDUCTION 


Exit plane data were reduced using the FORTRAN program 
plane.f. RVC3D output the following non-dimensional values 


for each grid point in the file fort.3: 


q,= 
re 
U 
eee 
* Poo 
Vv 
oe cal 
| ae 
pw 
Ci 
Poo 
ews a. 2 2 Ww 
fa | es =|— —— —. 
Polec.? ~4\ C Gane os ke: 





Computational output was manipulated in plane.f for 
comparison with experimental results. The following table 


provides a coordinate system comparison. 


Plot3ada 


Camvestan Cartesian 


Xie) ae Z(1) - q2 
ya) y(i) - q3 


zZ2(1) - w x(1) - a4 
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Transformation to cylindrical coordinates was conducted 
using the grid location information contained in the RVC3D 
fea anput file fort.1. 


6 (i) ee | cd) oP) 
(7 Ae 


All angles (8) were referenced to the y axis shown in the 
blade geometry plot of Figure 3. Flow angle computations were 
completed using the non-dimensional gq values described above 
as follows: 


Flow Angle ¥ (dimensional): 


I 


Va=-(u)sin@+(v)cos@ V,=(u)cos6+(v)sinO6 V,=w 


Zz 


v-tan{ | 
V. 


Flow Angle f_: 


_ PV,e __ | _ 
“ena pC, = come sinO+ (as) cos8 ee D.C, 





Direct comparison was possible with pressure ratio 
measurements. The ratio between local, and upstream reference 


total pressures was computed as follows: 





2 ae 


Co aa oh 
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Mach number waS obtained using the isentropic 


relationship between Mach number and temperature ratio. 











or T al 2, 2eee 
C= eet +Q, + 
I" TTy-ay TF 3g,(% q3+qa) 
T_y(y-1) 1 poe 
ae ds ~ Q2 +G3 +7 
ae qs 5 aq, 2 3 a) 

sm 
M= (72-4 } 2 

AG 7 


The program plane.f is included here for reference. It 
is an updated version of the program written by Thomas [Ref. 
4). The new version includes a calculation of flow angles 
after conversion to cylindrical coordinates for proper 


comparison with experimental results. 


CREKEKEKKEKEEEEKREKREKEKEEKEKRKEEREKKKEKERREKEKRKEREREEKEEKREREREKREEKEKEREKEKKKEKE 


plane. f reads rvc3d files & writes ascii files for plotting 
exit plane mach number am flow angles 
(comparison with experimental results) 


qaaQqgaaana 


unit 1 = input xyz file 

unit 3 = input q file 

unit 4 = output residual history 

unit 7 = output (k #1) comter, theta, pt/p0, Mach, Angle 
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e: unit 8 = output (k #2) counter, theta, pt/p0, Mach, Angle 
SG unit 9 = output (k #3) counter, theta, pt/p0, Mach, Angle 
c unit 10 = output (k #4) counter, theta, pt/p0, Mach, Angle 
e unit 11 = output (k #5) counter, theta, pt/p0, Mach, Angle 
E 


CER KEKKEKEKEKKEREREEREREERERKEREERRREKRERKERERRKKEEKREKRKRKEREREREKEKERKKKE 
parameter (ni=150,nj=31,nk=65) 
integer um 
real x(ni,nj,nk) ,y(ni,nj,nk) ,Z(ni,nj,nk), gama 
real qq(5,ni,nj,nk,resd(5000,5) ,ql,q2,q3,q4,q5,m(5) 
dimension kk(5) ,ang(5) ,theta(5) ,tt0(5) ,vtnd(5) ,vrmd (5), 
1 vznd(5) ,ptp0(5) ,deg(5) 


k-values are hard-wired below (10%, 25%, 50%, 75%, 90%) 
(use span.f to obtain %-k equivalency) 


Q20Q0 


data kk/21,27,33,39,45/ 
CEEKKEKKEKEKKEKKREKREEKEKEKEKRKEKKEKEKREREKEKRRRRRRREKRKREKREKKRKKEKREKRKEEKKEKKKKEK 
c read g71id coordinates 
CEEKEKKKEKEKEKKRKKEKEKEKRKEKREKRKEKREKEEKERKREKRERKEKEKEEKREKRKREREKEEREREKRRREKEER 

read (1,*)im,jm,km 

read (1,*) (((x(i,3,k) ,i=1, im) ,j=1, 3m) ,k=1,km) , 

1 (((y(1,3,) ,1=1, 1m) ,J=1,Jm) ,k=1,kn), 

2 (((Z(1,3,K) ,1=1,1m) ,J=1,)m) ,k=1,kn) 
CEBERKERRREKKEEKKEKKEEKEKEKKKEEKEKKKEKKKKEKKKKEEEKKEKKKKKKKEKKERERKEKEEEKEKEKKKKEEEE 
e: read restart file 
CEKKKKKEKEKRKEKEKEKREKERKKKRERKEKKKEREERKEEKKEERREKKEKEKKEREEREKEREKRKEKEKKKKKEER 

read (3,*) imax, Jmax, kmax 
read (3,*) fsmach,alpha,re, time 


c 
icheck=1abs (im-imax) +iabs (jm-jmax) +iabs (km—kmax) 
if (icheck.ne.0) then 
write (6,610) im, jm, kn, imax, jmax, kmax 
stop 
endif 
S 
read (3,*) ((((qq(1,1,),k) ,i=1,im) ,j=1,jm) ,k=1,km) ,1=1,5) 
bs 


6 additional residual data 
read (3,*)1tl,ii1,phdeg,gama,om,nres, dum, dum,dum,dum 
read(3,*) ((resd(nr,1) ,nr=1,nres) ,1=1,5) 
CEKKKEKKEKRKKEKKEEKKEKERERKEKREKEKEKEKKEREKEKEKKEREREEKKKERERREKEEEKKEKEKEKKKKKKKKE 
S Pressure side of exit wake cut... 
CEEKEKEKRKEKRKEKE KEKE KEKRKEKKEKEKEEKKKKERREEKKEKKEREKEKKEKREKEEEKKKEKKKKEKEKEKKKKEKE 
k=kk (3) 
‘=1 
do 10 j=jmax,1,-1 
um=um--1 
do 15 1=1,5 
k=kk/1) 
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Q 


Q 


Q 


Q 


©) 19) 19) (0) 1@, 


Rename desired q values for use 


qi=qq(1,1,3,k) 
q2=qq (2,1, ,k) 
q3=qq(3,1,3,k) 
q4=aq(4,1,),k) 
q5=qq(5,1,5,k) 


Campute pt/p0O ratio 
ptpo (1) =gama* (gama~1) *q5+ (gama-gama**2/2) /ql* (q2**2+q3**2+q4**2) 
Compute T/Tt ratio and mach number 


tt0 (1) =gama* (gama-1) /ql* (q5-0.5* (q2**2+q3**2+q4**2) /q1) 
m(1)=((1/tt0(1)-1) *2.0/(gama-1) ) **0.5 


Compute Theta and Non-D velocities theta, r ard z 


theta (1)=asin(y(i,3,k)/(2(1,3,k) **2+y (1,3 ,k) **2) **0.5) 
vtnd (1) =-q4*sin (theta (1) )+q3*cos (theta (1) ) 

vrnd (1) =q4*cos (theta (1) )+q3*sin (theta (1) ) 

vznd (1) =q2 

deg (1)=theta (1) *57.296 


Compute flow angle (degrees) using vtnd and vznd 


ang (1) =atan2 (vtnd (1) ,vznd(1))*57.296 
15 continue 
write *“7,300)um,deg(1) ,ptp0(1) ,m(1) ,ang(1) 
write (8,300) um,deg(2) ,ptp0(2) ,m(2) ,ang(2) 
write(9,300) um, deg (3) ,ptp0(3) ,m(3) ,ang(3) 
write(10,300)um,deg(4) ,ptp0(4) ,m(4) ,ang(4) 
write(11,300) um, deg (5) ,ptp0(5) ,m(5) ,ang(5) 
10 continue 


Repeat above on other side of wake cut at imax 
jJ=1 not used (duplicate location across wake) 
jJmax not repeated (same location as j=31 at i=1) 


k=kk(1) 
ql=qq(1,1,3,k) 
q2=aq(2,1,3,k) 
q3=q4(3,1,3,k) 
q4=aq(4,1,3,k) 
qo=qi4(5,1,),k) 
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ptpo (1) =gama* (gama-1) *q5+ (gama-gama**2/2) /ql* (q2**2+q3**2+q4**2) 
tt0 (1) =gama* (gama-1) /q1* (q5-0. 5* (q2**2+q3**2+q4**2) /q1) 
m(1)=((1/tt0(1)-1) *2/ (gama-1) ) **0.5 
theta(1)=asin(y(i,j,k)/(Z(1,3,k) **2+y (i,j,k) **2) **0.5) 
vtnd (1)=-q4*sin (theta (1) )+q3*cos (theta (1) ) 
vind (1)=q4*cos (theta (1) )+q3*sin (theta (1) ) 
vznd (1) =q2 
deg (1) =theta (1) *57.296 
ang (1)=atan2 (vtnd(1) ,vznd(1) ) *57.296 
25 continue 
write(7,300)um,deg(1) ,ptp0(1) ,m(1) ,ang(1) 
write (8,300) um,deg(2) ,ptp0(2) ,m(2) ,ang(2) 
write (9,300) um,deg(3) ,ptp0(3) ,m(3) ,ang(3) 
write (10,300) um,deg(4) ,ptp0(4) ,m(4) ,ang(4) 
write (11,300) um,deg(5) ,ptp0(5) ,m(5) ,ang(5) 
30 continue 
CEE KK KKK KKKKKKK KEKE KKK EKKEKRKKKK EKER EKER KK EKER KKKKKKKKKKKKKKKKS 
e residual history output to unit 4 
CEKKKKKEKKKKKEKKKKKKEKKKKKKKKKK KKK KKKKKKKKKKKEKKKKKKKKKEKKKKKKKEKKKKKKKKE 
write(4,310) 1, (resd(1,1),1=1,5) 
do 40 j=2,nres 
it=10* (j-1) 
40 write(4,310)it, (resd(j,1) ,1=1,5) 
CER KEK KEKKKKEKKKKKKKKKKKKKKKEKKKKKKKKKKKKK KKK KK KEK KEKE KKK KKKKKKKKKKKKKKE 
300 format(i5,5f8.3) 
310 format(i5,5(1x,e10.3) ) 
610 format (’ kkkkk warning KKKKKS | / 
aI ‘ im, jm, kn, read from input’,3i5,’ do not match’,/, 
2 ‘ im, jm, kn, read from restart file’ ,3i5) 
stop 
end 
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